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Abstract 

Within a transport code of BUU type the production of (p mesons in the reactions 
Ni + Ni at 1.93 A-GeV and Ru + Ru at 1.69 A-GeV is studied. New elementary 
(yT) [ reaction channels pN(A) and 7rN(1520) — > (j)N are included. In spite of a 

substantial increase of the (j) multiplicities by these channels the results stay below 
the tentative numbers extracted from experimental data. 
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1 Introduction 



Meson production in heavy-ion collisions is considered a sensitive probe of 
the reaction dynamics [1]. Particularly interesting is the production of heavy 
mesons at such a kinetic energy which is insufficient to create the respective 
meson in a free, single nucleon-nucleon collision. Secondary and cooperative 
processes are needed to accumulate the energy for the creation of the heavy 
meson. Therefore, the study of the subthreshold particle production is thought 
to reveal many aspects of the behaviour of strongly interacting matter includ- 
ing the equation of state [2]. 

Focusing on heavy-ion collisions at bombarding energies in the range of one 
A-GeV the production of strangeness degrees of freedom attracted much at- 
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tention, both experimentally and theoretically. For instance, the recent mea- 
surements seem to comply with a noticeable mass splitting between the 
and mesons [3] , as predicted by various effective (chiral) models for pseu- 
doscalar mesons [4]. 

In contrast to the comparatively rich body of experimental data on open 
strangeness meson production, the information on mesons with hidden strange- 
ness, i.e. the 0, is scarce. Only a few measurements [5,6] for the reactions Ni 
-I- Ni at 1.93 A-GeV and Ru -I- Ru at 1.69 A-GeV have been reported. In these 
measurements a strongly limited amount of the phase space was accessible. 
Extrapolations to the full phase space point to a surprisingly large production 
cross section. The results of the analysis [7] for the reaction Ni + Ni [5] let one 
recognise that the calculations based on a transport model underestimate the 
data by more than a factor of ten. (Notice that refined re-analyses [6,8] of this 
reaction has shifted the previously reported yield [5] to a larger value.) In the 
calculations [7] the mesons are solely produced in the elementary processes 
BB 0BB and ttB — > 0B, where the symbol B denotes the baryons N and 
A. Very recently the analysis of the production within another transport 
model has been repeated [9] with the same production channels confirming 
the previous finding [7] . The discrepancy of these calculations with the data 
[6,8,10] has inspired us to study whether further elementary channels might 
essentially contribute to the (j) meson production. 

In Ref. [11] the role of three-nucleon collisions for (p production was investi- 
gated. It turned out that these processes are well described by a sequence of 
two two-body processes inside the dense medium. However, to obtain the cor- 
rect yield, one has to include all intermediate particles which can be created 
inside the medium. This concerns especially the p mesons which have a rather 
large production cross section and a low threshold energy when colliding with 
baryons. 

Furthermore, even at subthreshold bombarding energies of about 1.5 up to 2 
GeV per nucleon, up to 30% of the nucleons arc excited to A and N* reso- 
nances which have a lower threshold for meson production than ground state 
nucleons when colliding with other particles. It is the aim of this study to in- 
vestigate the effects of these additional production channels which have not 
been studied so far. To be specific, we consider here the elementary reactions 
pB ^ B0 and 7rN(1520) ^ N(/). 

The study of the meson production is also of interest concerning the K~ 
production. Due to the strong decay channel (f) — > K'^K~ a substantial feed- 
ing of the K~ channel could not be ruled out. This is particularly tempting 
as the analysis [6] points to a substantial fraction of stemming from (f) 
decays. In the near future also the decay channel 4> ~^ e+e~ becomes exper- 
imentally accessible with the HADES detector at the heavy-ion synchrotron 
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SIS at GSI/Darmstadt [12]. Due to the large ss quark content, mesons are 
expected to interact much weaker with the non-strange baryonic matter in 
accordance with the Okubo-Zweig-Iizuka rule (cf. [13]). Nevertheless, notice- 
able in-medium effects, e.g. mass shift and broadening, are expected at higher 
temperature and density [14,15]. 

Our paper is organised as follows. In section 2 the cross sections of the above 
mentioned elementary hadron reactions are estimated within a one-boson ex- 
change model. Sections 3 and 4 describe the resulting multiplicities when 
implementing these new reaction channels in a transport code. In section 5 we 
compare with available experimental results. The summary and the discussion 
can be found in section 6. 



2 Elementary cross sections 

We are going to describe heavy-ion collisions within the Boltzmann-Uehling- 
Uhlenbcck (BUU) formalism. At beam energies around 2 GcV per nucleon 
many resonances are generated. Therefore, there are many possibilities to pro- 
duce mesons by two-body collisions. In Refs. [7,9] the cross sections of the 
collisions of the most abundant particles, NN — > NN </>, NA — > NN0, AA — > 
NN (f), ttN N0, and ttA N0, have been used. Since the cross sections 
are not well known a one-boson exchange model has been employed to es- 
timate them [7,16]. The parameters of the corresponding Lagrangians have 
been fixed by decay widths and cross sections measured at different energies. 
Selected new data near threshold became recently available [17]. As mentioned 
above, employing these cross sections within transport codes one gets (j) meson 
production rates which seriously underestimate the extrapolations made from 
the measurements [5,6]. 

Channels with p mesons as incoming particles have not been included in 
Refs. [7,9]. Since the meson decays into the non-strange sector overwhelm- 
ingly via the irp channel (with 13% branching ratio, see [18]), both mesons are 
expected to play an important role for the (f) production. Other channels are 
substantially weaker. 

2.1 4> "meson production in p-baryon collisions 

Let us consider first the production in p induced reactions on the nucleon 
and the A, i.e. pN (pN and pA — > (pN. Figure la shows the correspond- 
ing tree level diagram hitherto not yet considered. We calculate the cross 
sections within the one-boson exchange model with the effective interaction 
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Fig. 1. Diagrammatic representation of processes contributing to the reactions (a) 
p (N,A) ^ N, (b) TT N* ^ N and (c) np cf). 

Lagrangians 



■'ttATA ■ 



ttNA 



m 



ipTiP^ ■ d^ir + h.c, 



(1) 
(2) 
(3) 



Here ^0 is the nucleon field, is the Rarita-Schwinger field for the A reso- 
nance, and 71, pI, are the meson fields while m^r and stand for the pion 
and p meson masses, respectively. The symbol f denotes the isospin Pauli 
matrix, and T is the isospin 1/2 — > 3/2 transition matrix. The diagrams in 
Fig. 1 are to be understood as guidance to construct Lorentz invariant ampli- 
tudes from effective interaction Lagrangians. They must be supplemented by 
additional elements to achieve a reasonable description of elementary hadron 
processes, for which data are available. In this line of arguments we intro- 
duce multipole form factors for off-shell pions with four- momentum at each 
vertex which, e.g. for the nNN vertex reads 

y v^ttnn) ~P-k J 

with n — 1, while for p mesons we use a corresponding expression with n = 2. 
The cross section for the pA — * (pN reaction is regularised by the Peierls 
method [19]. We adopted the coupling constants fi and cutoff parameters Aj 
for the 71 NN and ttNA vertices from the Model 11 in Table B.l of Ref. [20]. 
For f^^p^ and Ajf^^ we used the values given in [7] . 



4 





n 




A~ 


AO 


A+ 


A++ 


p 





2 








1 


3 




1 


1 





2 


2 







2 





3 


1 









Table 1 

Isospin factors for the various channels of the reactions pN <^N and pA 
The isospin averaged differential cross sections can be written as 
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with s being the square of the center-of-mass energy and Pin(out),CM the incom- 
ing (outgoing) three-momentum in the center-of-mass system. The quantity 
I Alp represents the squared matrix elements summed over the final and av- 
eraged over the initial spin and isospin states. These quantities are given by 
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for the pN channel and 
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Kt>N\ 




X (pa • Pn + mAmN) 



'Pa-PnY 2 
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{{Pp ■ P<l>f - mlml) (7) 



for the pA channel. Here pn, Pa, Pp, and denote the four- momenta of 
the respective particles. In the case of the pN channel, pn and p'^ refer to 
the four-momenta of the incoming and outgoing nucleon, respectively. The 
factors Ai are given by ApN->,iyN = {UNNf7rpcp/m^rn^)KNNiPl)Kp^iPl) and 
ApA^4>N = {f7rNAfnp,t>/m^rn^)F:;^NA{pl)Kp4>iPl)- The total cross sections are 
obtained by numerically integrating Eq. (5) over the sohd angle Q. 

Cross sections of the specific isospin channels can be deduced from the isospin 
averaged cross section by multiplying with an appropriate isospin factor. These 
isospin factors, which can be expressed in terms of Clebsch-Gordan coefficients, 
are given in Table 1. 



2.2 (f) formation from the np channel 



The cross section of the irp process has a very narrow resonance structure with 
a width of Ttot — A A3 MeV and a branching ratio of T^p/Ttot — 0.13 in the 
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7rp channel. The formation cross section, depicted in Fig. Ic, of the meson 
in the TT+p" or the 7r~p"'" channel is correspondingly 



a — 



2pi(Vi-m^)2 + r,V4 



which gives a maximum cross section of 9.5 mb at resonance peak position. 
The n^p'^ channel is isospin forbidden. 



2.3 Effects of higher resonances 



At bombarding energies of about 2 GeV per nucleon baryons can be highly 
excited. Consequently the production threshold is reduced by their excitation 
energy if mesons are created in collisions with other particles. To estimate 
the importance of such processes we test the influence of the collision of a 
TT meson with the N(1520) 3/2^ resonance. The corresponding diagram is 
exhibited in Fig. lb. This resonance is chosen because it decays with a large 
branching ratio of about (20 ±5)% [18,21] into the p channel, although this 
channel would be closed for the nominal masses of the p and N(1520). This 
hints to the fact that a large part of the production cross section proceeds via 
virtual p propagation as described by the diagram in Fig. lb. This assumption 
may not be applicable to higher resonances which can emit p mesons nearly 
on shell. 

There are several possibilities to describe the pNN(1520) vertex by using sim- 
ple effective interaction Lagrangians: 

jC.pNN{1520)^goi^Nfi>ir,2oP^ + h.C, (9) 
^pNN(1520)^—'li^NrY'il^l52o{di.Pu.-dp,p^) + h.c, (10) 

^pNN{1520) = ^'0JVT(7"^'0152O ^a^/iP/? + h.C. . (11) 

nip 

The second Lagrangian is related to the vector dominance approach (see e.g. 
[22]) while the last one follows from constituent-quark models [23]. We fix the 
coupling constants by assuming that each Lagrangian alone would provide 
the partial width to the p(two-pion) decay channel. Such a procedure supplies 
upper limits on the coupling constants. We obtain the values go — 7.5, /i = 
10.5 and /2 = 68, respectively. With these values the Lagrangians (9) and 
(10) give roughly the same cross section near threshold. Our coefficient /i 
is larger than the corresponding coefficient in Ref. [22], which is caused by 
the monopole form factor used there. In Fig. 2 we present the cross sections 
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Fig. 2. Production cross sections for cj) mesons for several binary reactions as a 
function of the relative kinetic centre-of-mass energy in the entrance channel. The 
cross sections are calculated for the nominal masses of the resonances. Especially 
cross sections for pA collisions having different masses can considerably deviate 
from the curve shown. 

obtained with the Lagrangians (10) and (11). The large coupling constant 
/2 needed for the Lagrangian (11) leads to a much larger cross section for 
the 4> production than the Lagrangian (10). This value of /2, however, differs 
from the value estimated from the quark model [23]. As our Lagrangians are 
nothing more than a covariant parametrisation of the interaction the coupling 
constants are generally energy dependent. The larger relevant energy in 
production compared to the N(1520) decay could be the reason for the different 
estimated values of the coupling constants. In the following we use the cross 
section from the Lagrangian (10) as standard in our BUU model calculations. 
We also performed calculations using the cross section calculated with the 
Lagrangian (11) in order to estimate the uncertainties of the meson yield 
from the tt + N(1520) channel, see Fig. 2. 
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2-4 Summary of cross sections 



Summarising the results of this section we show the calculated individual cross 
sections in Fig. 2 as a function of the kinetic energy of the relative motion 
in the entrance channel. Obviously, the p B channels are very important at 
low energies. We also studied the reaction tttt — > 07r which, however, is not 
important because of the large threshold energy which can rarely be overcome 
by the pions. For the (p production via the BB and ttB channels we use the 
same cross sections as found in Ref. [7]. 



3 BUU calculations 

To calculate meson production in a heavy-ion reaction we use the version of 
the BUU code introduced in Ref. [24]. The production of mesons proceeds 
via the reactions discussed in the previous section. To increase the statistics 
the perturbative method is used, i.e. if in a two-particle collision the threshold 
is overcome, the (j) meson is created and weighted with its production proba- 
bility. However the production process itself and its subsequent collisions do 
not influence the motion of the other participating particles. To simplify the 
calculations an isotropic cross section is assumed for each elementary hadron 
reaction and in the case of a three particle final state the energy distribution 
is taken from the phase space. 

In the transport model all nucleon and delta resonances are included up to 2.2 
GeV. Furthermore, A and E baryons as well as the mesons 7r,r), p,a,u!, 
and are explicitly propagated. The properties (masses, widths, branching 
ratios) of the baryon resonances have been fitted [25] to the one- and two- 
meson production data available from pion-nucleon reactions. The production 
cross sections for resonances in nucleon-nucleon collision are derived from data 
on meson production in nucleon-nucleon reactions. 

Furthermore, momentum dependent potentials for the miclcons arc included. 
The same potentials are also used for the A particles. As standard potential 
we employ such one which delivers a soft equation of state with an incom- 
pressibility parameter of 215 MeV. Less is known about in-medium properties 
of the mesons. While in Ref. [15] a very small mass shift has been found, 
we use for exploratory purposes the prediction by Hatsuda et al. [14] for the 
density dependence of the 4> meson mass: 

n 

ml = m<^ - 0.025m^— . (12) 
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Fig. 3. Number of produced p and (p mesons as a function of time for the reaction 
Ru on Ru at beam energy of 1.69 A-GeV and impact parameter h = 2.8 fm. The 
lower pair of curves is calculated by switching off the p production due to resonance 
decays. 



The second term is treated as the mean field of the meson. 

The interaction of the mesons with surrounding nucleons is relatively weak. 

The clastic cross section is estimated to be in the order of 0.5 mb, however, a 
strong absorption cross section is expected via the reaction channel 0B — *KA. 
Following Ref. [26] we use a value of 6mb ■ exp(— 3.3e/GeV) with e being the 
kinetic energy in the entrance channel. Similar values are employed in Ref. [7] . 

Using these potentials we have calculated the production cross sections for 
the reactions Ni on Ni at 1.93 A-GeV and Ru on Ru at 1.69 A-GeV bombarding 
energy. The results are reported in Table 2. 
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Table 2 

Multiplicities of 4> mesons per central event in Ni ± Ni and Ru ± Ru reactions. 
The first 5 lines give the results from the special channels indicated. The symbol B 
comprises the nucleon and the A particle. Lines 7 and 8 are the yields observed in 
the two detector subsystems. 

4 Influence of model parameters 

In the first five lines of Table 2 the contributions of the different production 
channels are hsted. We notice a substantial increase by a factor of 2.5 of the cf) 
multiplicities when including the new p channels. The tt + N(1520) reaction 
channel calculated with the Lagrangian (10) does not essentially contribute 
to the production. However, when using the cross sections obtained from 
the Lagrangian (11) the (p meson yield from the tt + N(1520) channel rises 
by more than an order of magnitude to 6.6-10-^ (7.7-10"^) for the reaction Ni 
+ Ni (Ru + Ru). This indicates that the tt ± N(1520) channel might not be 
negligible for the production, but to clarify this requires a better knowledge 
of the pNN(1520) vertex. 

In order to test the dependence of the (p meson yield on the equation of state we 
carried out also calculations with different potentials. For instance, employing 
a hard momentum dependent potential with an incompressibility parameter 
of 380 MeV only two third of the yield of the soft equation of state has 
been obtained for the Ru on Ru reaction. The reduction of the (p multiplicity 
concerns especially the BB and ttB channels while the pB and irp channels 
are less influenced. Obviously, to overcome the larger potentials the collision 
partners use up a large amount of kinetic energy which then lacks for the (p 
production. Similar effects are known for kaon production, see Ref. [3]. We 
have also found that momentum independent potentials enhance the number 
of the produced (p mesons. 
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Fig. 4. Relative number of produced (j) mesons as a function of the density at creation 
point for the reaction Ru + Ru at 1.69 A-GeV beam energy 

The time dependence of the number of produced p and mesons is exhibited 
in Fig. 3. The time origin is chosen such that at 5 fm/c the nuclei touch. 
The maximum p number is obtained at complete overlap (13 fm/c). However, 
the maximum density is attained very early at 9 fm/c. At this time the 
production reaches its maximum value and ceases at 15 fm/c. 

In the calculations it turned out that indeed most of the p mesons are gen- 
erated by resonance decays; a substantially smaller part stems from the tttt 
annihilation process. To demonstrate the importance of the baryon resonance 
decay into p mesons for meson production we tentatively switched off these 
decay channels. Then the p number is diminished to 25%. This in turn reduces 
drastically the number of produced 0's by a factor of two. 

To estimate the importance of the in-medium effects of the mesons in dense 
nuclear matter it is interesting to see at which densities the 0's are created. 
In Fig. 4 we display the number of mesons as a function of the local density 
at the creation point. The production reaches its maximum at around 2.5 
times normal nuclear density while the average density is about twice normal 
nuclear matter density. Nevertheless the effect of the mean field in Eq. (12) is 
not dramatic. The yield increases by about 10% only. The effect is stronger 
in the BB channels. Furthermore, the slope parameters of the momentum 
distribution (see discussion in the following section) arc reduced by about 5 
MeV. One should be aware that the in-medium effect of mesons has to be 
treated together with the effects of the K mesons. If e.g. in addition the sum 
of the K~ and K+ masses is reduced by 200 MeV (see [4]) the decay time of 
the meson is considerably shortened and could be smaller than 5 fm/c [15]. 
Such an effect would decrease the detection probability due to the occurrence 
of kaon rescattering. This effect is not treated here. Due to the 0B— >KA cross 
section about 30% of the created mesons are reabsorbed. 
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In Ref. [7] a multiplicity of 3-10~^ has been obtained within a transport 
model treatment for Ni on Ni collisions. This value is smaller than ours, even 
if we do not count the contributions from the p and N(1520) channels. Besides 
the use of different mean fields, one has to observe that in our calculations 
the effect of the increasing decay widths due to the dropping K~ mass has not 
been included. The formation from the K'^K" channel is negligible as was 
already shown in Ref. [7]. 



5 Comparison with experiment 

The FOPl collaboration has measured meson production in the reactions 
Ni + Ni at 1.93 A-GeV and Ru + Ru at 1.69 A-GeV [5,6] via the decay 
channel — > K^K~ . The detector system allows to identify the original 0's 
in two separated parts of the momentum space. The analysis [5] based on 
data of identified pairs within the acceptance of the central drift chamber 
(CDC) at target rapidity was recently supplemented [6,10] by taking into 
account also the HELITRON forward Plastic Wall subsystem of the FOPl 
detector which covers a rather limited phase space near midrapidity. A full 
phase-space extrapolation of the observed K='= pairs from the very limited 
acceptance regions of both detector systems can be performed only by relying 
on assumptions on the spectral shape of the meson distribution. Using the 
information from both detector systems, CDC and HELITRON, and assuming 
a thermal distribution of the mesons, dN/dp oc exp {—E/Tes), one has been 
able to estimate the total yields given in the last fine of Table 2. Thereby 
a rather low temperature of Teflf — 59 (72) MeV has been obtained for the 
reaction Ni + Ni (Ru + Ru). 

Our calculations have been carried out for central events comprising about 9% 
(15%) of the total cross section in the symmetric collisions of Ni (Ru) nuclei. 
Despite the gain of mesons due to the large contribution of the p mesons 
the calculated numbers are still smaller than the </> multiplicities extrapolated 
from experiment. For completeness we list also in Table 2 the yields within 
the acceptance regions of the CDC and HELITRON subdetectors. 

In our calculations we find that the momentum spectrum of the (f) mesons is 
not compatible with a thermal distribution. The distribution is stretched out 
in longitudinal direction. We can attribute two different slope parameters T± 
and Tj| for the perpendicular and longitudinal distribution, respectively. Fit- 
ting independently the distributions near midrapidity and at small transverse 
momenta we find for the Ru + Ru (Ni + Ni) collision T± = 90 (100) MeV and 
3][ = 113 (120) MeV, respectively. The rescattering of mesons increases the 
longitudinal temperature by about 10 MeV while the transverse temperature 
is less affected. 
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Fig. 5. Predicted momentum distributions of (p mesons in the centre-of-mass system 
for central Ni + Ni collisions for two centre-of-mass angles 0. The error bars on 
the symbols reflect the statistical uncertainty in our Monte-Carlo calculation. The 
smooth curves show fits to a thermal model with temperature values of 120 MeV 
for Gem = 30" and 100 MeV for Gem = 90'', respectively. 

To illustrate the effect of the anisotropy of the meson emission we exhibit in 
fig. 5 the momentum distribution at two different centre-of-mass angles. Both 
distributions are fitted by a thermal model with angular dependent tempera- 
ture parameters. 

In this respect it is interesting to observe that the different channels provide 
very different angular distributions. The mesons which are created in BB 
collisions have a nearly isotropic angular distribution. However, in the ttB 
and pB channels the resulting angular distributions are significantly forward- 
backward peaked. This anisotropy caused by the meson channels determines 
the total momentum distribution. This leads to the fact that the yield in the 
backward lying CDC detector is larger than that in the HELITRON detector 
which is sensitive at midrapidity. This finding seems to be in contradiction to 
the preliminary analysis by the FOPI group [10]. Also results by Aichelin and 
Hartnack [9] within a quantum molecular dynamic transport code support this 
behaviour of the different channels. As in their calculations the BB channels 
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Fig. 6. Transverse momentum spectra of electron-positron pairs from 0— >e~''e~ 
decays at three different centre-of-mass rapidities y shown by full symbols. Corre- 
sponding open symbols indicate those electron-positron pairs which can be observed 
within the geometry of the HADES detector. The efficiency of the electron and 
positron detection is not accounted for. The lines are drawn to guide the eye. 

play a more important role they obtain a nearly thermalised distribution 
though the cross section is essentially smaller there. 

The anisotropy can be studied by dedicated measurements with the HADES 
detector. To get an idea of the expected spectral shapes we display in Fig. 6 
the transverse momentum distributions of the electron-positron pairs which 
stem from meson decays for three rapidity slices. The open symbols in Fig. 6 
indicate the pairs which enter the HADES detector^. One observes that the 
filtered e"'"e~ spectra resemble much the original (j) — > e"'"e~ spectra depicted 
by the full symbols in Fig. 6. This suggests that the original distribution 
can be recovered by the detector. 



We thank R. Kotte for providing us with the routines of the HADES acceptance. 
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6 Summciry and discussion 



In summary we have presented an analysis of (f) meson production in central 
heavy-ion collisions at beam energies of 1.69 and 1.93 A-GeV. As new channels 
we have included pB dN, Trp — and vr N* ^ 0N and copious p production 
via the resonance decays B* pN. Despite of enlarging the multiplicities, as 
follows within a transport code of BUU type, the experimentally found values 
[6,10] are still underestimated. This conclusion is not affected by the theoreti- 
cal uncertainties regarding the tt -|- N(1520) channel. A momentum dependent 
soft mean field for baryons and a weak mean field for the mesons have been 
employed. We have not considered a possible change of the decay width of 
the meson as the momentum dependence of the in-medium properties are 
poorly known at present. 

With respect to the potential importance of the meson for the K~ channel 
a dedicated measurement with enlarged statistics would be highly welcome 
to resolve the puzzle. It should be emphasised that the meson becomes 
now also accessible via the e+e" decay channel in HADES measurements [12]. 
The branching ratio between the electromagnetic and the hadronic decay will 
be sensitive to the in-medium properties of mesons. Thus, the K^K^ decay 
channel should be addressed in future experiments with both the up-graded 
FOPl apparatus and the HADES detector due to the precise tracking system. 
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